STATISTICAL DERIVATION OF BASIC 
EQUATIONS OF DIFFUSIONAL 
KINETICS IN ALLOYS WITH 
APPLICATION TO THE DESCRIPTION 
OF DIFFUSION OF CARBON IN 

AUSTENITE 

V. G. Vaks, I. A. Zhuravlev 

National Research Center "Kurchatov Institute" , 123182 Moscow, Russia 

Moscow Institute of Physics and Technology (State University), 
117303 Moscow, Russia 



Basic equations of diffusional kinetics in alloys are statistically derived using the 
master equation approach. To describe diffusional transformations in substitution 
alloys, we derive the "quasi-equilibrium" kinetic equation which generalizes its ear- 
lier versions by taking into account possible "interaction renormalization" effects. 
For the interstitial alloys Me-X, we derive the explicit expression for the diffusivity 
D of an interstitial atom X which notably differs from those used in previous phe- 
nomenological treatments. This microscopic expression for D is applied to describe 
the diffusion of carbon in austenite basing on some simple models of carbon-carbon 
interaction. The results obtained enable us to make certain conclusions about the 
real form of these interactions, and about the scale of the "transition state entropy" 
for diffusion of carbon in austenite. 



1 INTRODUCTION 

The problem of development of an adequate theoretical description of diffusion in al- 
loys attracts interest from both fundamental and applied standpoints, see, e. g. [T]-[T5]. 
Presently, this description is usually based on the phenomenological theory of diffusion 
in multicomponent systems developed by Onsager many years ago [6] . Phenomenological 
kinetic coefficients of this theory are calculated using various simplified models with pa- 
rameters estimated empirically [TJ-[I]. However, these empirical models have usually no 
consistent theoretical justification, and their relation to interatomic interactions, as well 
as possibilities of applications to other alloy systems, are typically not clear. 

One of important problems in this field is the strong concentration dependence of the 
diffusivity D of carbon in austenite [TJ-[5]. This dependence causes complications in the 
kinetic analysis of various diffusion-controlled reactions in steels [I] , and several empirical 
models have been suggested to describe this dependence [TJ-[1]. However, possibilities to 
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employ these models for predictions of D at temperatures T < 1000 K (where many 
important phase transformations occur but D can not be directly measured as austenite 
is unstable here) or under partial substitution of iron by other metals are, generally, not 
clear [5]. 

The consistent statistical description of the diffusional kinetics in alloys can be based 
on the master equation approach [7]- [T5]. This approach enables us to express the phe- 
nomenological kinetic coefficients, such as the mobility M a of an a-species atom, via in- 
teratomic interactions in an alloy These interactions can be estimated using either some 
microscopic models [16]- [18] or ab-initio methods [TJ1 [19], [20] . As the level of accuracy 
and reliability of ab-initio calculations is steadily increasing, this microscopic approach 
seems to be prospective for non-empirical calculations of diffusivitiy. 

At the same time, previous considerations of diffusional kinetics in alloys based on the 
master equation approach were usually restricted by discussions of only simplest models 
and approximations or some particular problems [Zl El O [HI HBj- On the other hand, 
few more general discussions [TUl [T5] included many not necessary complications and 
restrictions which can hinder the understanding of the results. Therefore, the first aim 
of this paper is to present a clear and general derivation of basic equations of diffusion 
in alloys based on the master equation approach, for both substitution and interstitial 
alloys. 

To this end, in Sec. 2 we first discuss the vacancy-mediated kinetics under diffusional 
transformations in substitution alloys. This problem has been considered in Ref. [15j . 
and some equations derived in Sec. 2 have been already presented there. However, the 
full derivation of these equations has not been given in [15], while the similar derivation 
in [10] includes a number of complications and inaccuracies. In Sec. 2 we also discuss the 
methods of computer simulations based on the microscopic equations proposed, including 
some generalizations of the earlier-discussed "equivalence theorem" [TUl E3 which greatly 
simplifies such simulations. In Sec. 3 we consider the interstitial alloys Me-X and derive 
the general statstical expression for the diffusivity D of an interstitial atom X which has 
a simple and physically transparent form. This expression includes only microscopic pa- 
rameters which can be estimated using either theoretical models or ab-initio calculations. 
We also generalize this microscopic expression for D for the case of multicomponent alloys 
(MeiMe2...)-X with several species atoms in the metal sublattice. 

The second aim of this work is to apply the results of Sec. 3 to microscopically 
treat the above-mentioned problem of diffusion of carbon in austenite. This treatment 
described in Sec. 4 is based on the microscopic model of carbon-carbon (C-C) interactions 
in austenite suggested by Blanter [TB] which supposes a strong "chemical" repulsion at 
short C-C distances Rcc an d a purely deformational (or "strain-induced") interaction at 
longer distances -Rcc- We show that some natural generalizations of this model enable us 
to describe both thermodynamic and diffusional characteristics of carbon in austenite at 
the same level of accuracy as that achieved in phenomenological models [U [2j [3j HI [36] . At 
the same time, the microscopic approach used enables us to make a number of conclusions 
about the type of carbon-carbon interactions and about some physical features of diffusion 
of carbon in austenite. The main results of this work are summarized in Sec. 5. 
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2 EQUATIONS OF VACANCY-MEDIATED KINET- 
ICS FOR DIFFUSIONAL TRANSFORMATIONS 
IN SUBSTITUTION ALLOYS 

2.1 General equations for mean occupations of lattice sites 

First we present the necessary relations from Ref. [15] with some extensions and com- 
ments. We consider a substitutional alloy with (m+1) components p' that includes atoms 
of m different species p = pi,p2, ■ ■ - Pm and vacancies v : p' = {p, v} . The distributions of 
atoms over the lattice sites i are described by the different occupation number sets {nf } 
where the operator is 1 when the site i is occupied by a p'-species component and 
otherwise. For each i these operators obey the identity ^ , n\ = 1, so only m of them 
are independent. It is convenient to mark the independent operators with Greek letters 
p or a : (nf )i n d ep = rf, while the rest operator denoted as is expressed via nf : 

nf=(l-E4 C 1 ) 

p 

Note that both and are the projection operators: 

(n*) 2 = nt, nf< = 0, « = ( 2 ) 

For dilute alloys, it is convenient to put u h" in (pQ) to be the host component, e.g., /i = Fe 
for the dilute BCC Fe-Cu-u alloys discussed in [UJ [15] and used below for illustrations. 

In terms of all operators nf the total configurational Hamiltonian H l (for simplicity 
supposed to be pairwise) can be written as 

Ht = \ E v ^ «■ ( 3 ) 

After elimination of the operators according to ([I]), this Hamiltonian takes the form 
H* = E + <P P < + H int , H mt = <°< n ° ( 4 ) 

pi pa,i>j 

which includes only independent while constants E , (p p , and "configurational inter- 
actions" are linearly expressed in terms of the couplings V? q in §3$), in particular: 

v pa = fypa _ yph _ yha + yhh^ (5) 

The fundamental master equation for the probability P of finding an occupation number 
set {<} = £ is [12]: 

dP(OM = Yl^ w ^ ^Piv) - w(v, 0^(0] = ^ (6) 
»? 

where W(£,rj) is the 77 — ^ transition probability per unit time. If we adopt for prob- 
abilities W the conventional "transition state" model [H], we can express the transfer 
matrix S in ([H]) in terms of the probability of an elementary inter-site exchange ("jump") 
pi ^ vj between neighboring sites % and j: 

Wff = <n^JJ exp[-/3(^. - E% vj )] (7) 
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where (5 = 1/T is the reciprocal temperature, E^£ V j is the saddle point energy, -Ep? ■ is 
the initial (before the jump) configurational energy of a jumping atom p and a vacancy, 
and the pre-exponential factor icyj can be written as 

u$f = u in e X p{ASg vj ). (8) 

Here u pv is the attempt frequency which, generally, is expected to have the order of 
magnitude of a mean frequency of vibrations of a jumping atom in an alloy, and AjS^ - is 
the entropy difference between the saddle-point and the initial alloy states. This difference 
is mainly due the difference of atomic vibrations in the saddle-point state (supposed to be 
locally equilibrium so that such thermodynamic notion as the entropy can be applied to 
it) and the initial state. At high temperatures T under considerations (actually, already at 
T > 9/2 where Or, is the Debye temperature [22]), this entropy difference can be described 
by the classical expression: 

AS s pv p = 31n(</<) (9) 

where u p n and are certain average frequencies of vibrations of a jumping atom in the 
initial and in the saddle point states, respectively, see, e. g., [2H]- As frequencies ui^ in the 
saddle-point configuration can notably soften with respect to u p n , the entropy difference 
AS p p can be expected to take significant positive values. For example, for the Fe-Cu-f 
alloys with the Debye frequency cu^f ~ 6 ■ 10 13 sec -1 [23], Soisson and Fu (SF) [T4] found: 
u>t{ ~ 80 u%, u e c {{ ~ 30 u$. It implies: AS§f v ~ 4.5, AS* P V ~ 3.5, cog ~ w £/4, 
Wq u ~ w^ e /3, though so high values of AS SP and u D /u S|) can be somewhat overestimated 
due to inaccuracies of estimates |14j . 

The saddle point energy E^ v j, generally, depends on the atomic configuration near 
the ij bond. We describe this dependence by the model of SF [T4] supposing the saddle- 
point energy to depend only on occupations of lattice sites / nearest to the center of bond 
ij (denoted by l l J n ): 

= E < n i= E i + K- ( 10 ) 

. jij 

y 3 t — >>nn 

Here E v h is the saddle point energy for the pure host metal, while the operator A?- de- 
scribes changes in this energy due to a possible presence of minority atoms near the 
bond: 

where z h nn is the total number of nearest lattice sites I for each bond (being z b nn = 6 
for the BCC lattice), A p p = {e p — e p h ), while (e p and e^) are the microscopic parameters 
of pairwise interactions calculated by SF using ab-initio methods [14]. Note that our 
definitions of A p - and A p p differ by sign from those used in [T4j and |15j . 

The interaction parameters V? q , e v h and A^ in Eqs. (jSJ) and ffTT]) can be calculated 
by ab-initio methods. For the Fe-Cu-f alloys it is illustrated by SF [T4]. Theoretical 
calculations of factors ic'tjj in Eq. (J^J) are more difficult due to the presence of the 
entropic factor exp(AS' 5 ) in Eq. (JHJ). However, values of uj e JJ can be estimated from 
experimental data about the self-diffusion and the diffusion of isolated atoms in a host 
metal, as described by SF [T4] . 

As n\ in Eqs. ([TJ, (j2J), and §6§ are the projection operators obeying Eqs. (j2J), the most 
general expression for the probability P = P{nf} in ([6]) can be written in the form of 
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the generalized Gibbs distribution []IJ [HI [12] 

P{<} = exp |>(n + J] W - #int - ^nt)] , (12) 

pi 

= ? E K<™° + lH Kk< n >i + ■■■ ( 13 ) 



Here parameters Af (which are both time- and space-dependent, in general) can be called 
"site chemical potentials" for p-species atoms; they are related to local chemical potentials 
pf and p^ of p-species and host atoms as Af = (pf — pj 1 ) [21] . H int in (fT2j) is the same 
as in dl]); the parameters in (Tl3|) (also depending on both time and space) describe 
possible renormalizations of interactions; and Q is determined by normalization. 

As discussed in detail in [12], for the usual conditions of phase transformations cor- 
responding to the absence of external fluxes of particles or energy (that is, when the 
alloy is a "closed" but not an "open" statistical system), the effects of renormalizations 
of interactions can be expected to be insignificant, thus one can put in (fl2l) : 



hint = 0. (14) 

There are at least two reasons to expect the validity of Eq. (|14p for the transformations 
in closed systems. First, this relation holds both before and after the transformation. For 
example, it is true before an initially equilibrated alloy is quenched from a higher tem- 
perature Th to the lower temperature TJ corresponding to another equilibrium phase (or 
phases), and it is also true after the new equilibrium state at T = T\ is reached. Therefore, 
there is no driving force drawing the correlation parameters h p i '"° in the distribution (I12p 
away from their Gibbs values h%"'J = 0. Second, the parameters h?'"J in the distribution 
(JT2J describe mainly the short-range order. After a change of external conditions, such as 
temperature, this short-range order is established relatively fast, for a time of the order of 
one interatomic exchange time r a , while the time for completing microstructural evolution 
under phase transformation is usually much longer: t 3> r a [SI [TUl [15] • Therefore, possible 
fluctuative violations of relation ([141 at small t < r a are not important for the whole 
evolution. 

When relation (IT4"|) is obeyed, Eq. (fl2l) takes the form: 

P{< } = exp [f3(Q + K< - Him)] (15) 

pi 

which will be called the "quasi-equilibrium Gibbs distribution" . 

Note, however, that for the essentially "open" systems, such as alloys under irradia- 
tion [24, 25J or an alloy with an external atomic flux imposed [JTJ, the quasi-equilibrium 
relation f[T4l can be significantly violated. Important cases of such violations can be the 
phase transformations accompanied by significant fluxes of vacancies into the transforma- 
tion region, for example, the precipitation in Fe-Cu-f alloys where these fluxes arise due 
to the strong trapping of vacancies by the copper-based precipitates [HE] . In such cases, 
the effective vacancy-copper interactions can notably vary with the evolution time, and 
the presence of significant renormalizations h^ u ^ can be expected. Possible methods 
of calculations of these renormalizations have been discussed in [TO], [Til H3] • 

Multiplying Eq. ([6j) by operators and summing over all configurations {n?"}, we 
obtain the set of equations for mean occupations of sites ("local concentrations") cf = 

«>: 

dcfjdt = (n>S) (16) 
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where ((...)) = X^n^C-O^i^J} m eans averaging over the distribution P, for example: 

= K> = ZXP{nJ}. (17) 

For what follows it is convenient to mark the minority atoms by Greek indices a, /3, ... 
Then index p in Eqs. (T7|)- (|T2|) is a or h, corresponding to a minority or a host atom, while 
index p in Eqs. ( FT2|) -( [T6|) is a or v, corresponding to a minority atom or a vacancy. Then 
the system of equations (|T6|) can be explicitly written as follows: 

dc a r i 

= E E w « ex p - /^sy -{<-»• #J ex p[^( fi - ^//)] 

«} i(0 

dc v r i 

= E E E < f K n i ex P O 3 ^* - ^4') - {* i)J expL9(" - ^//)] (18) 

«} p j(0 

where symbol j(i) means summation over sites j being nearest neighbors of site i, and 
p is h or a, while H e ff=H e ff{n^} is the effective Hamiltonian for statistical averaging of 
expressions in square brackets: 

Heff = - e k< + \ E < = Kr + Ki)- (19) 

p,fc pa,kl 

Here t>^ can be considered as the full effective interaction, which for simplicity is supposed 
to be pairwise, just as the "true" interaction in the hamiltonian fll]). The operator 
E™ v j in ([18]) (describing that part of the initial configurational energy which depends on 
occupations of sites % and j) can be expressed in terms of formal variational derivatives 
of the hamiltonian ()3]) over n P i and n v i, H l pi = 8H 1 /5n P i and Hp i v j = 5 2 H t /Sn P iSn V j: 

^pi,vj = n pi-H p i + n V jH V j — n P in V jH piv - (20) 

where the third term corresponds to the substraction of the "double-counted" interaction 
between an atom p at site % and a vacancy v at site j. 

The main idea of further manipulations (analogous to those made in [101 [TTJ) is to 
reduce the averages of complex operators in square brackets in (fl8~|) to some more simple 
averages having a clear physical meaning. To this end, in sums over all configurations 
{n^} in (I18p we first perform summation over all possible occupations of only two sites, 
% and j, belonging to the ij bond under consideration. Due to the presence in (1T51) of 
the projection operator nfn?-, the summation reduces to putting n\ = 5 qp and = 5 pv in 
the nf- and n^-dependent exponential factor exp Y multiplied by this projection operator, 
where 

Y = /3(E™ vj -H eff ). (21) 

Note that the saddle-point energy E pi p v j, according to its definition (jTUJ), does not contain 
operators of occupations of site i and j. Therefore, the common factor u e JJ exp(/3f2 — f3E p ^ v -) 
will be skipped for brevity in Eqs. (I23|) - (l26|) . 

For what follows it is convenient to restore formally the summation over all occupation 
number sets in ([TBI , including all values of and nj. To this end we can introduce 
the operator n^n 1 ^ in the summand. Since this projection operator is nonzero only when 
all nf and rf- are zero, the summation with this factor over all possible occupations of 
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sites i and j is equivalent to omitting all nf- and nj-dependent terms in the exponential 
exp Y. Therefore, the result of the summation can be written as 

J2 <n) Y = n 1 nl J ex P ( Y pi + Y -J + Y p^ + Y ) ( 22 ) 

O rr O rr 



% ' 3 t ' 3 



Here Ypi^ Y v j cincl V j are the variational derivatives of the operator Y over the relevant 
occupation numbers: Y pi = 5Y/Sn^, etc. The first term of the exponential in Eq. (j22p 
corresponds to the contribution to the sum (1221) of the term in Y linear in nf but not in 
rij, the second, to that of the term in Y linear in nj but not in nf, and the third, to that 
of the term in Y linear in both vR and n"". 

Let us first consider in (|22p the term with p = a and express all operators nf in the 
expression (TSUI) for E™ ■ via the independent nf using Eq. (jTJ . After the above-mentioned 
putting in this expression nf = 1, = 1, the exponent Y in (12~T|) takes the form 

F = [ E far - < + E fa^ p - *s» fc ) < + E fa* + - - H *ff] ■ (23) 

Using Eqs. (Tl9i) and ([5]) we can explicitly write the relation (122]) at p = a as follows: 
<n v j exp Y = ^2 nfnj exp {/? [j>5 + «J)nf + £ (V + V,f ) - - Kg 



n p ,n a . n p ,n a Ip 



/p 

where the quantity 

<i = {Vt ~ VH h ) (25) 

can be called "the kinetic interaction" for a p-species atom (as it influences only mo- 
bilities but not thermodynamic properties [U]). Note that the vacancy concentration 
c\ = (n?) in real substitution alloys is very small, hence all nf can be neglected in the 
statistical averages which enter Eqs. (Tl8l) . Therefore, terms with nf in (T2^1) correspond 
actually to the minority atoms with p — (3 ^ v. In writing Eq. ([24]) we also used the 
above-mentioned considerations that for the usual conditions of phase transformations, 
the significant renormalizations h\f can be expected only for the vacancy-atom interac- 
tions, while for the interactions between different atoms, the analogous renormalizations 
hff are not essential. 

In the case p = h, Eq. (122]) is simplified as the operator Y in (12~T]) depends only on 
the independent operators n^ but not on . Therefore, terms Yu and i/w )t y in (122]) are 
absent, and the exponent is reduced to (Y v j + Y). Making the same manipulations as 
described above, we obtain for this case instead of 



<n) e W Y=J2 <n) exp {/? [$>2 + «>f + £ faf* + V^f) - V™ 

n p ,n tT n p ,n tJ 1/3 I 

+ (*;-£*#»?)-#«//)]}■ (2 6 ) 

ip 

Substituting relations (1241) and (126]) in ([18]) we can express the derivatives dc^/dt 
via some statistical averages. In writing these averages we can take into account that 
the interaction renormalizations h int in (j!2p are present only for the vacancy-atom terms 
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Kj n l n j which include the vacancy occupation operators and thus can be neglected. 
Therefore, the quasi-equilibrium distribution ( Fl5|) can be used in calculations of these 
averages. To simplify formulas, in treatment of interaction renormalization terms h^nf 
in Eqs. and ( )26|) we will use the mean-field approximation (MFA), replacing each 
operator nf in these averages by its mean value cf. Therefore, each term h v ^n^ is 
replaced by h^f cf , which corresponds to replacing the vacancy site chemical potential 
by its "renormalized" value A^: 

(k - e *y»r) k?<f ) ■ (27) 

Z/3 Z/3 

In Sec. 2.2 we show that at cf values close to unity, when the interaction renormalization 
effects can be expected to be most significant, the approximation fl27|) becomes exact. 
Hence, one can expect it to be sufficiently accurate at all cf. We can also note that a 
similar MFA treatment of interaction renormalization effects have been used by Nastar 
et al. [TT| [T5] , and comparing their MFA and kinetic Monte Carlo results, these authors 
found the MFA accuracy to be sufficient for treatments of the renormalization effects. 
Using relations (T2^|) - (f27]) we can write Eqs. (fl8i) in the concise form 



dct/dt = ^7^(0#-£X) 



m 

dcydt = Y,[^{l h ,b1 j + Y,^ b ^)-{^J}] (28) 

m 

which will be called the "quasi-equilibrium" kinetic equations (QKE). The term r y pv in 
these equations (where p is a or h, i. e. denotes a minority or host atom) is the effective 
exchange rate p f^v for a pure host metal. This term can be written in the form similar 
to Eq. ©: 

7p , = < / exp(-/3^) (29) 

where u e JJ is the same as in (jBJ, while E%! is the effective activation energy which is 
expressed via the saddle point energies E\ in ffTTl) and interactions V? q and hf? in 



and (TT9j) as follows: 

E h a v c = e£ - J2(v£ h + V$ h ) + V£ (31) 

j 

where nn means "nearest neighbors". Comparing these expressions to the analogous 
activation energies E p ^ cMC used in the kinetic Monte Carlo approach [3] and given by 
Eq. (2.5) in pj, we find: 

jpav _ rpav _!_»?«». jphv _ rphv (on\ 

^ac — ^aCjMC ^ U nni ^ac ~ ^a^MC V°*J 

where v av is the same as in (TT91 . The difference between and E^ MC arises because 
in the statistically averaged QKE f[28l) . the transition probability (J7|) is averaged over the 
distribution (fT2"l) . For the inter-site exchange ad ^ vj, it leads to an extra Gibbs factor 
exp(— f3v° v ) with from ffT9l) in the averaged probability. 
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In this connection we note that in their study of diffusion in dilute Fe-Cu alloys 
with a thorough taking into account the correlation effects for this vacancy-mediated 
diffusion [27], SF obtained for an effective activation energy of a copper atom the value 
(E^ v ) eff ~0.47 eV [251125]. This is very close to the value E^ v =0AA eV which follows 
from Eq. (I3"U1) at hf™ = after substitution of couplings V? q used by SF. It may imply 
that for the diffusion in dilute Fe-Cu alloys, our statistical averaging with the quasi- 
equilibrium Gibbs distribution ([15]) can rather accurately describe the correlation effects 
mentioned [27], while the interaction renormalization effects for this diffusion are not very 
significant. 

The quantities in (125]) (to be called "correlators") are certain averages of site 
occupations which describe influence of minority atoms in vicinity of the bond ij on the 
pi^vj jump probability: 

^. = (^nJexp[^/3« + ^)<- E ^ A >"]> ( 33 ) 

al n i—jij 

ct, i — i"nn 



where A£, is the same as in (fTTl) . and is the same as in 

Finally, the quantities and rj" in (1251) can be called the "site thermodynamic 
activities" for vacancies and a -species atoms, respectively, because they are related to 
the site chemical potentials Af in (1151) and the renormalized site chemical potential A^ in 
fl27J) as 

C = exp(/3A7); r£ = exp (/?A?) , (34) 

that is, similarly to the relations between conventional thermodynamic activities and 
chemical potentials. 



2.2 Calculations of site chemical potentials and correlators 

To find explicit expressions for site chemical potentials Af = Af(cj) determined by Eqs. 
ffTT]) . and for correlators 6^- = t^j(ck) determined by Eqs. (1331) . we should use some 
approximate method of statistical physics, such as the MFA or cluster methods [T2] . 
As discussed in detail in [151 EH EH], employing the MFA for calculations of chemical 
potentials Af in real alloys often leads to great errors, while the pair cluster approximation 
(PCA) usually combines simplicity of calculations with a high accuracy, particularly for 
dilute alloys. As an illustration (used also below for interstitial alloys), we present the 
PCA expressions for Af in a binary A-B-v alloy with host atoms h=A and minority atoms 
a=B (skipping index B at Af = Aj and cf = Ci for brevity): 



(35) 
(36) 



Here the function or is expressed via the Mayer function = [exp (—(3vij) — 1] 
or f?- = [exp (— (3vff) — 1] for the potential i>y = u?. B or t>^ B defined in (jSJ) as follows: 



9ij 
9h 



R 



'.i 



2fi/[R ij + l + Mc i -c j )) 
-- {[1 + (cj + c^] 2 - Aacjfijifij + 1)} 



1/2 



(37) 
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For a multicomponent alloy A-B!-...B m -f , the PCA methods of calculations of site chem- 
ical potentials Af are described in [2T] . 

In calculations of correlators (1331) . we first consider for simplicity the case of 
configuration-independent saddle-point energies when differences K p a and in Eqs. 



TOD, (TIID and (133 



are zero and correlators 



by are independent of the kind of a 



jumping atom p. Using Eqs. (jSJ) and the identity 



exp (xnf) = 1 + nff{x), f(x) = {f - 1), 



(38) 



which follows from we can rewrite Eq. ( I3"3"j) as 



fc=0 h^...l k OL X ...OL k 



(39) 



where we set 

= + /?„«) (40) 

with /(a;) from (1381) . while fcf in ([39]) is the total number of sites with nonzero values of 
potentials (u^ + For example, for the nearest-neighbor or next-to-nearest-neighbor 
interaction models in a BCC lattice [T5] . we have fc t = 14 or fc 4 = 20. 

In finding averages in f[3"9"j) we should consider that the functions /" in Eqs. (I3"9~j) 
and (1401) for real alloys are typically rather large. For example, for the BCC Fe-Cu-v 
alloys considered in [15] we have: f(flui) ~ 5 and f{fiu2) ~ 1 (where interaction 
Mi and Ui corresponds to the nearest and next-to nearest neighbors). Thus the main 
contributions to sum (1391) come from averages of products of many different operators 



nf corresponding to well-separated and weakly correlated sites I. In particular, for the 
BCC lattice, these products (even for the nearest-neighbor interaction model) include 
terms with the neighbors from first to tenth, most often third and fourth. Correlations 
of occupations of so distant sites should typically be weak. Therefore, using the simple 
MFA that neglects such correlations, in calculations of averages (1391 should, generally, be 
adequate, unlike calculations of chemical potentials Aj mentioned above. 

In the MFA, each operator nf in Eq. (139]) is replaced by its average value cf . Hence 
the correlator bij can be explicitly written as 



c^Cj exp 



[X>(i+£/r< 

1=1 a 



1=1 



(41) 



When the differences A^, in Eqs. (1111) and (1331) are nonzero, the correlator in Eq. 
(13"3"|) can be calculated by the same way as b^ in (1391) - (14T1) . The difference arises only for 
sites I = V£ n adjacent to the ij bond, for which the factor /" defined by Eq. (T4T1 is 
replaced by an analogous factor f££ defined as 



rap 

h 



/A 



(42) 



where A^ 



A^ ji, and 5 l ^ is unity when I = l l J n and zero when I ^ l^ n . Therefore, 



the correlator is given by Eq. (l4~Tj) with replacing each by / ( 



ap 
IA '■ 



'.I 



cfcjexp 



hi 



ap a 
IA C 1 



1=1 



(43) 
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Finally, let us make a remark about the MFA-type approximation (|27p used above in 
the derivation of QKE (1281) . If we don't use this approximation, Eqs. ( |28|) . instead of the 
correlators (|43|) . include similar correlators 6f- differing from by the presence of 
additional interaction renormalization terms Ph^nf in the exponents: 

kt 

% = c h lC ) exp J2 ^ { 1 + E [ ex P + - A "/ - " X ] c " } ( 44 ) 

i=l a 

where we write the function f(x) from (138|) explicitly. For the diffusional transformations 
under consideration, the interaction renormalization effects seem to be most significant in 
those space regions where the local concentration cf is close to unity. It is illustrated by 
the case of precipitation in Fe-Cu alloys where these effects arise due to the strong trapping 
of vacancies by the Cu-based precipitates for which cf u is close to unity [IS]. For such cf , 
the argument of logarithm in (j44p is reduced to the single exponent, and the relation fl27|) 
becomes exact. Let us also note that for the simulations of diffusional transformations 
based on Eq. (1251) . the details of vacancy distributions are actually insignificant due to 
the "adiabaticity principle" and the "rescaling of time" procedure discussed in Sec. 2.3. 
Therefore, the approximation (|27|) appears to be sufficient for using in such simulations. 



2.3 Reducing kinetic equations (1281) to those for some direct 
exchange model 



The QKE (1281) can be used for modeling of most different phase transformations, in 
particular, of processes of precipitation which attract great attention in connection with 
numerous applications [3U1 ED E2J • However, in the original form (j2"81) these equations 
are not suitable for using in computer simulations due to very small values of vacancy 
concentration c v in real alloys. As atomic exchanges pi ±= vj take place only with a 
vacancy, this smallness leads to a great difference in the relaxation time r between atoms 
a and vacancies v: r a ~ r v /c v ^> t v . It is illustrated by the presence of vacancy activities 
= exp(/3A^) in the right-hand side of the QKE (T2SJ. This activity is proportional to 
the vacancy concentration c^ which is a general relation of thermodynamics of dilute 
solutions illustrated by Eq. (136|) . Therefore, the time derivatives of mean occupations are 
proportional to the local vacancy concentration, c^ or cj. It is natural for the vacancy- 
mediated kinetics and leads to the strong inequality between r Q and r v mentioned above. 
Therefore, the type of temporal evolution for atoms and vacancies is quite different, which 
makes the direct numerical solving of Eqs. (l2"gj) for c?(t) and cf (i) to be unsuitable and 
time consuming. 

At the same time, the inequality r Q 3> t v enables us to use the "adiabatic" approach 
used in many fields of physics, including the well-known Born-Oppenheimer approach 
in the quantum-mechanical description of motion of atoms in molecules and solids. In 
this approach, the effective driving force for a slow motion is obtained by its averaging 
over a rapid motion. Therefore, to fully describe the slow motion, only few averaged 
characteristics of the rapid motion are needed. In the quantum mechanics, this is the 
appropriate electronic energy ("electronic term") calculated at the fixed positions R» of 
atoms. In our problem it means that at the given atomic distribution {cf}, the local 
vacancy concentration c v t adiabatically fast (i. e., for a time r v ~ c v r a r a ) reaches its 
"quasi-equilibrium" value c^{cf } for which the right-hand side of the second equation 
(128]) vanishes. Therefore, discarding small corrections of the relative order c" C 1, we 
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can approximate this equation by its adiabatic version: 



(45) 



which can be called "the adiabaticity equation" for the vacancy activity Solving this 
linear equation for we can, in principle, express it via c". Then substitution of these 
Ci(Cj) into the first Eq. (12 8p yields the QKE for some equivalent direct-atomic-exchange 
(DAE) model. 

To illustrate these considerations, we first consider the models with the configuration- 
independent saddle-point energies. For such models, the parameters A? in ( ITT]) are zero, 
the correlators = 6^ do not depend on the kind p of a jumping atom, and the 
adiabaticity equation ( I45p takes the simple form: 



Ihv 



+ £7«tf) .7'}' 



(46) 



If we let 1/vi denote the first term in the square brackets (14"6"|) , then the difference in these 
brackets takes the form (u^ 1 — vj 1 )- Hence, a solution of Eqs. ( l4"6j) is given by z/j being a 
constant independent of the site number i (although possibly depending on time as well 
as on temperature and other external parameters): 



Vi = H I (y hv + lavVl 



(47) 



Relation ( 14 7p determines the above "quasi-equilibrium" vacancy distribution c^{c"} which 
adiabatically fast follows the atomic distribution {cf}. Substituting it in the first Eq. 
([28]) we obtain an explicit kinetic equation for atomic distributions {cf} for which the 
evolution of the vacancy distribution is characterized by a single parameter u(t) being a 
"spatially self-averaged" quantity: 



dcf/dt = Y b v v (t) l™lhv (Vj ~ V?) 

m 



(4f 



Equations (I48|) can also be rewritten in the form used for DAE models [12J: 

dcf/dt = Y, M t 2 smh[/3(AJ - \f)/2] 

m 

+ M ? 2 sinh[/3(A° + Af - Xf - Af )/2] 



(49) 



where the generalized mobilities M^ 9 , which describe the inter-site exchanges a^h and 
are given by 



M. 



ah 
ij 



l*vlhvv{t) %exp[/3(A? + A")/2] 



(50) 
(51) 



M?f = lavl p v v{t) % exp [f3(\f + X] + Af + Af )/2]. 

Comparing these expressions with the expression (32) in [TU] which describes the mobility 
Mff in an alloy with the nearest-neighbor direct-exchange rate lij—l Pq , we see that Eqs. 
(j50p and (jSTl) correspond to a DAE model with the effective direct exchange rates 



lah = Icxvlhv V{t)\ lap = lavlfr f (*)• 



eff 



(52) 
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As u(t) in ( H7|) is proportional to c?, the effective DAE rates (|52p are by a factor c„ smaller 
than the vacancy exchange rates j pv , in accordance with the above-discussed adiabaticity 
relations. 

For more realistic models with configuration- dependent saddle-point energies, the ba- 
sic adiabaticity equation (I45p for vacancy activities Q can not be solved analytically 
in general, and hence either numerical or some approximate analytic methods should be 
used. For the first-principle model of Fe-Cu-f alloys developed in [13], such approxi- 
mate treatments made in [15] have shown that the equivalence relations fj48|) - fl52|) usually 
preserve their form, but the correlator by is replaced by some other quantity, bf^ or bj? 
Physically, the possibility to reduce the vacancy-mediated kinetics to the equivalent direct 
atomic exchange kinetics is connected with the above-mentioned fact that in the course 
of evolution of an alloy, the distribution of vacancies adiabatically fast follows that of the 
main components. Therefore, it can be assumed that this equivalence is actually a general 
feature of the vacancy-mediated kinetics, while for more general models, the correlators 
bij in (T5T)|) can be replaced by some other expressions with similar properties. 

The function v(t) in Eq. (152 p determines the rescaling of time between the initial 
vacancy-mediated exchange model and the equivalent DAE model fj48|) - fl52|) . Temporal 
evolution of this DAE model is actually described by the dimensionless "reduced time" 
t r related to the real time t by the differential or integral relations 



where r^ h = l/j^k nas ^ ne meaning of the mean time of an atomic exchange a h, 
while t r has the meaning of an effective number of such atomic exchanges. 

The form of the function t(t r ) in ( 15 3 p depends on the boundary conditions for vacan- 
cies adopted in simulations. In particular, if we adopt the "vacancy conservation" model 
for which the interaction renormalization effects can be expected to be insignificant, we 
can use Eqs. (j3Z|) and (T36|) to express the local vacancy concentration c\ via v{t) and 
cf(t r ). Then the vacancy conservation condition: £\ c^(t,t r ) — N v — const, can be used 
to explicitly find the dependence t(t r ). However, taking into account a possible creation 
of vacancies at various lattice defects (grain boundaries, dislocations, etc) suggested in 
the kinetic Monte Carlo (KMC) simulations [T4"l 126] appears to be more realistic. Then 
the dependences t(t r ) can be found from comparison of results of the DAE-based simu- 
lations described above to the appropriate KMC results, as illustrated in [15], and these 
dependences seem to be rather simple and universal. In more detail, applications of Eqs. 
(Hgp - (T53"P to studies of precipitation in concrete alloys will be described elsewhere. 



3 EQUATIONS FOR DIFFUSION OF INTERSTI- 
TIAL ATOMS IN INTERSTITIAL ALLOYS 



In binary interstitial alloys Me-X where X is an interstitial atom, in particular, in iron- 
carbon steels, diffusion of atoms X is realized via thermo-activated jumps of these atoms 
between their interstitial sites ("pores"). Therefore, this diffusion can be described by the 
general equations of Sec. 2 for a particular case of a substitution binary alloy X-t> which 
consists of atoms a=X and vacancies v in the crystal lattice of pores, with the "host" 






(53) 
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atoms h being vacancies v. The total configurational Hamiltonian here includes X-X 
interactions between atoms X, but not X-t> and v-v interactions. Therefore, only V xx 
terms are nonzero in formulas (JHJ-flSJ): 

V?« = V™ + 0, V t f = V t f = 0. (54) 

The only meaningful index a=X is usually skipped below, for example: cf = q, V^ x = 
thus the effective Hamiltonian fjTTJ]) takes the form: 

H eff = ~ ^ K^i + H int , H int = g v ii n i n r ( 55 ) 

i ij 

The mean occupation q = (n,) of pore i by an atom X is related to the local chemical 
concentration by the relation depending on the geometry of pores [2J, e. g. c=x/(l — x) 
for a uniform austenite structure MeX c with the FCC lattice of octo-pores. 

An important principal feature differing the diffusional kinetics in the interstitial Me- 
X (i. e., substitution X-u) alloys from that in the substituton A-B-v alloys is the validity 
for Me-X alloys of relation (|14|) . that is, the absence of interaction renormalization effects. 
It follows, first, from the physical considerations presented after Eqs. (|14p and ( 1T5|) and, 
second, from the thorough analysis of interaction renormalization effects for A-B-f alloys 
made by Nastar et al [HJ [13]. These authors found these effects to be described by the 
terms antisymmetric with respect to interchanging A and B atoms, (h AB — hfj A ), which 
vanish in a binary X.-V alloy where A=B=X. Therefore, the diffusional kinetics in Me-X 
alloys can be described by the quasi-equilibrium relations (fT4"]) and (fT5]l . 

For a uniform Me-X alloy, the site chemical potential Aj=A in fl55]) coincides with 
the thermodynamic chemical potential fix, unlike substitution alloys where the analogous 
quantity A Q , as mentioned in Sec. 2.1, is equal to the difference (/i Q — /i^). To show it, 
we generalize Eqs. (21)-(24) and (40)-(43) of Ref. [2T] to the case of interstitial alloys 
Me-X. The quasi-equilibrium Gibbs distribution (fT5"|) and the generalized grand canonical 
potential Q g for the effective Hamiltonian (155]) have the form 

P = exp [(3(Q g + ^2 ^aiiiai ~ H int ) , (56) 

ai 

n g = -Tin ^ exp [fff KiUai - ^mtj , (57) 

{n ai } ai 

while the mean occupation q is related to {l g {\ a i} by the formula obtained by the 
differentiation of equality ( 157|) : 

Ci = (th) = -dn g /d\. (58) 
Therefore, if we define the generalized free energy F by the equality 

F = n g + J2^, (59) 

i 

then site chemical potential A« is related to F by the relations generalizing a similar 
relation for a uniform alloy: 

= dF/dd. (60) 

To relate Aj and Q g in Eqs. fl56p - (|6TI]) to the thermodynamic chemical potentials, we 
consider the case of a uniform alloy Me-X when and Aj in Eqs. (JHS1- (JSOD are independent 
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of i: Q = c, Aj=A. For definiteness, we discuss the austenite structure for which the total 
number of interstitial sites (octo-pores) is equal to the total number Nyi e of Me atoms. 
Then instead of the total thermodynamic potentials, Q g and F, it is convenient to 
consider the analogous quantities per one atom Me, Q and /: 

{l = tt g /N Me , f = F/N Me = n + Ac, c = N x /N Me . (61) 

Here Nx is the total number of atoms X, thus c is the mean occupation of an interstitial 
site, and, according to Eq. ([SO]) : 

\ = df/dc. (62) 

The quantities Q and A in Eqs. f l6T|) and ( 16"2"|) are simply related to the partial chemical 
potentials, /xx and /iM e , defined by the thermodynamic relations: 

fix = dF/dN x , // Me = dF/dN Me . (63) 

Substituting relations (I6"T|) for c and F = N Me f(c) in Eqs. ( 16"3"|) and taking into account 
Eq. fl62|) . we obtain: 

A = /X X , ^ = /iMe- (64) 

Thus, the quantity A or f2 in Eqs. flBTj) and (162]) has the meaning of the chemical 
potential of atoms X or atoms Me, respectively. 

Kinetic equation describing diffusion of atoms X in an interstitial alloy Me-X (treated 
as a binary substitution alloy X-v) can be obtained by substitution of relations (j54"|) and 
([55D into the first Eq. (EE}: 

dd/dt = ^2'j ij bf j [exp((3\ j ) - exp(/3Aj)]. (65) 

m 

Here we took into account that for a host atom h = v , the site chemical potential A^ 1 
in ( I2"gj) should be formally put zero, as illustrated by Eqs. (12"6"]) and (I2"8"]) . The jump 
probability 7^ and the correlator in fl65|l are defined by the relations analogous to 
(|29|) and (|33|): 

7jJ =^ / exp(-/3^), (66) 
fig = ((1 - n,)(l - % ) exp ( - /3 ^ Ajn fe ) }. (67) 

The pre-exponent u;^- in f !6B]) is determined by Eq. OH]) with replacing p — kX, while the 
activation energy E 1 ^. is reduced to the term E v h = E^ in ffTUl) . unlike the more complex 
expression (I3TI) in a substitutional alloy. Index "ij" at the quantities cu^J and E l J c in 
(I66p allows for a possible nonuniformity of an alloy; for a uniform alloy this index can 
be omitted. The quantity in (IBTI) is the analogue of Af- in Eqs. (JTDJ) and ffTTl) : it 
describes the change of the saddle point energy E§f v j for an inter-site X-atom jump i—tj 
due to the presence of another atom X at site k. 

The kinetic equation (165]) can also be written in a form analogous to Eq. fT4"9]) : 

da/dt = 2M ij sinh [/3(Aj - Ai)/2] (68) 

where the generalized mobility M^-, according to (165]) . is determined by the relation: 

= H3 h l ex P [/5(A< + \)I2\. (69) 
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In usual diffusion problems, the space dependence of functions q = c(rj), Aj = A(cj), 
and bfj = 6 x (q,Cj) in Eqs. (I65i) - (l69l) is supposed to be smooth. Therefore, variations of 
these functions under replacing q — > Cj (or r.j— )-rj= + where = rj — is the 
inter-pore distance) are small. Then the kinetic equation (165]) or (1551) can be expanded in 
powers of r^Vc. It yields the continuous version of this kinetic equation having the form 
of the "continuity equation" for the flow j of atoms X: 

dc/dt + divj = 0, J a = -J2 D ^{c)V u c . (70) 

V 

Here a and v are Cartesian indices, and the diffusivity D av is determined by the following 
expression: 

D av {c) = T au b x da x /dc, (71) 

r- = ^$>^ (72) 

m 

where ax = ax{c,T) = exp[/3X(c,T)] is the thermodynamic activity of X atoms, and 
bx = bx(c,T) is the correlator bfj in Eq. ( IBTj) at q = Cj = c. 

The site chemical potential A; in a binary alloy can be written as 

A, = \f + Af* (73) 

where A* d = Tln[cj/(1 — Cj)] corresponds to the ideal solution, and Af 1 * describes inter- 
action effects, see, e. g., ( 135]) . Therefore, for the ideal solution for which both Af 1 * in (173]) 
and Af in (167]) are zero, we have: ax = c/(l + c), 6x = (1 — c) 2 , bxdax/dc = 1, and 
Eq. (|7Up takes the form of a simple linear diffusion equation: 

dc/dt = J2K d ^Lc (74) 

with the concentration-independent diffusivity D^ v equal to in (172]) . However, when 
X-X interactions and are significant, the kinetic equation ([70]) is nonlinear, and 
the diffusivity D in (ITT]) should vary with the local concentration c = c x (r). 

For a uniform cubic alloy, such as austenite, tensor T au is reduced to a scalar 5 av ^a 2 
where 7 is given by Eq. (166]) (with omission of index ij) and a is the FCC iron (7-iron) 
lattice constant, hence D au = S au D. Then the diffusivity D, according to (ITT]) , can be 
written as 

D = ^a 2 bxdax/dc (75) 
where 7 is defined by Eqs. (166]) and ([8]) with the appropriate change of indices: 

7 = u e x ff exp(-f3E ac ), co e x ff = wxp (AS| P ). (76) 

For a uniform alloy with c, = c, the PCA expression (135]) for the local chemical potential 
A(c) is simplified [29] : 

A(c,T) = T^2;„ln(l - g n c), g n = 2f n /(R n + 1 + 2cf n ), 

n 

f n = exp (-Pv n ) - 1, # n = [1 + 4c(l - c)^] 1 / 2 (77) 

where z n is the coordination number, and v n is the configurational interaction for the nth 
coordination sphere. In the case of weak interaction, (3v n 1, Eq. ( TTT]) is reduced to the 
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MFA expression: A MFA = (Yl n z n v n) c. However, for the realistic values of interactions 
v n , such as those presented in Table 1 below, using the MFA can lead to significant errors 

Ill- 
Correlator bx{c, T) for a uniform alloy, according to Eqs. ([57)1 and (jUJ), can be written 

as 

b x (c, T) = (1 - c) 2 exp [ ]T ln(l + /* c )] , /* = exp (-/3A n ) - 1 (78) 

n=l 

where is the coordination number, and A„ is the saddle-point interaction for the 
nth coordination sphere of the saddle point considered. If these interactions are weak: 
/3 A n < 1, Eq. (EH]) takes its MFA form: 

b x (c, T) = (1 - cf exp (~PcJ2 4 P ^n) (79) 

n=l 

However, for the realistic saddle-point interactions A n , such as those shown in Table 2 
below, using the MFA can lead to significant errors, just as for A in Eq. (IT?]) . 

The microscopic relation (I73|) can be compared to various phenomenological models 
for diffusivity It can provide, in particular, the statistical expression for the phe- 

nomenological mobility Mcva introduced by Agren in his discussion of diffusion of carbon 
in austenite [1]- [3]. Comparing Eq. (T75l) to the definition of Mcva by equation (9) in [2], 
we find: 

M CVa = ja 2 b x a c / [c(l - c)V rn T] (80) 

where V m is the volume per atom Me. In more detail, the microscopic and phenomeno- 
logical descriptions of diffusion of carbon in austenite are compared in Sec. 4. 

Let us now consider a multicomponent interstitial alloy (MeiMe2-..)-X with several 
species atoms p in the metal sublattice, such as an Fe-Mn-C alloy. The interstitial sites 
will be marked by indices i, j and k, while the sites in the metal sublattice, by indices I, 
m and n. The total configurational Hamiltonian can be written in the form generalizing 
Eq. ©: 

Ht = \ll Wi + + E VSntf + \H V i» q m ( 81 ) 

ij p,il pq,lm 

where we again skip index X for an interstitial atom putting nf = rii, V^ p = V£, 
V-j X = Vij. As above, we discuss diffusion of only interstitial but not metal atoms, and the 
presence of vacancies in the metal sublattice is neglected. Therefore, occupation operators 
for host metal h can be expressed via those for the minority metals a similarly to Eq. 
(PQ): = (1 — Yl a n f)- ^ ne e ff ec tive Hamiltonian for statistical averaging, instead of 
(fT9~|) . takes the form: 

Heff= -J2 X ^~J2 X >1 

i a,l 

+ \ E W»i + E V > iTl l + l E (82) 
ij a,il ce/3,lm 

where v% = (V^ — V^) , while is related to in (1HT1) similarly to Eq. ([5j). 

Equations describing diffusion of atoms X can be again derived using Eqs. (15])- (|2^|) 
with appropriate generalizations and simplifications. In particular, the first equation (I18p 
here takes the form 

d*/dt= E^^-^Kexp^-^Af )-{i^j}]exp\p(n-H eff )]. (83) 

{n k , n?} 
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Here the saddle-point energy £g p , instead of Eqs. ([TUP and (fTTj), is given by the expression 

4f = E h + J2 Aj?n fc + A >f (84) 

A; 

where A^ and A^ are analogues of Ajjf in (1671) . Using relations (l20i) and (|26i) . we again 
reduce Eq. (183]) to the form (1651) . However, the activation energy i? ac in (1BB1) and the 
correlator bfj in (IBTj) now are defined as 

E ac = Eh-^ v iii 
i 

bf 3 = ((1 - n«)(l - n 3 ) exp ( - Aj?n* - £ A>f) ), (85) 

and Eqs. f[8"3"j) and flHBI) now include the statistical averaging over various distributions of 
a-species atoms in the metal sublattice. 

4 CALCULATIONS OF DIFFUSIVITY AND AC- 
TIVITY OF CARBON IN AUSTENITE FOR SIM- 
PLE MODELS OF CARBON-CARBON INTER- 
ACTIONS 

To calculate the diffusivity D according to the microscopic expression (1751) . we should 
use some theoretical model of X-X interactions in an alloy, for both the configurational 
interactions v n in (J77J) which determine the chemical potential A, and for the saddle-point 
interactions A n in (ITS]) which determine the correlator &x- For substitution Fe-Cu alloys, 
such first-principle model for both v n and A„ has been developed by SF [H], and simu- 
lations of precipitation in Fe-Cu alloys based on this model confirmed its adequacy and 
reliability [2] . For interactions of carbon in austenite, reliable first-principle calculations 
are still absent due to the well-known difficulties of taking into account magnetic inter- 
actions in 7-iron [19]. However, some simplified model of configurational interactions v n 
in austenite has been suggested by Blanter [IB] , and his estimates of these interactions 
are presented in Table 1 as . Below we use this model and some its extensions to 
investigate the concentration and temperature dependences of diffusivity D which follow 
from the microscopic expression (I75|) . 

Blanter used the model of purely deformational configurational interactions with the 
nearest-neighbor Kanzaki forces for all constants except the first one. The nearest- 
neighbor constant vf (which can not be adequately described by the deformational model 
due to the strong "chemical" repulsion at short C-C distances) was treated as a free 
parameter which was estimated from the fit of the carbon activity in austenite with 
respect to graphite, ajf gr , calculated with these to the experimental values. The 
quantity clq 9T is related to the "configurational" activity a>c = exp (/5Ac), where Ac = A 
is the chemical potential of carbon in austenite discussed in Sec. 3, by the thermodynamic 
relation [IB] : 

a c - 9r = a c exp(/3AG c - £ ' r ) (86) 

where AG c _9r = AGq~ 9T (T) is the difference between the thermodynamic potentials per 
carbon atom in a pure 7-iron and in graphite. The fit to experimental a c _9r (c, T) obtained 
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Figure 1: Dependence of carbon activity a?(7 9r in austenite with respect to graphite on 
the carbon concentration xq = c/(l + c) for various temperatures T. Dots correspond to 
experimental data presented in [16J. Solid curves are calculated using the PCA expression 
(|77p for A with the interaction constants v n from Table 1. Dashed curves are calculated 
using the phenomenological description of C-C interactions employed by Agren [36J. 

Figure 2: Fe-C phase diagram. Dots correspond to experimental phase boundaries. Solid 
curves show ferrite-austenite phase boundaries calculated using Eq. (177]) with v n from 
Table 1. Dashed curves correspond to the phenomenological calculations by Agren [36J. 

with the use for both Monte Carlo [16j, and the PCA [29J calculations of ac(c, T) the v% 
values, and for AGjf 9r (T), some experimental estimates, seemed to be quite satisfactory. 
It may imply that the simple model of Blanter [16] can serve as a basis for realistic 
descriptions of C-C interactions in austenite. 

Table 1. Configurational interactions v n (in kelvin) of carbon atoms in austenite 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


2R n /a 


110 


200 


211 


220 


310 


222 


321 


400 


330 


411 


420 


R n / Ri 


1 


1.41 


1.73 


2 


2.24 


2.45 


2.65 


2.83 


3 


3 


3.16 




12 


6 


24 


12 


24 


8 


48 


6 


12 


24 


24 


v%, Blanter [E] 


1334 


1961 


-487 


46 


46 


267 


-23 


-139 


58 


-12 


-23 


v n , this work 


1400 


1180 


-322 


46 


46 


267 


-23 


-139 


58 


-12 


-23 



In the present work, the configurational interactions v n have been estimated using a 
similar approach. However, in the fit to experimental values a? ( 7 9r {c ) T) we used for the 
function AG]7 9r in (1861) the interpolation of experimental data suggested by Agren [36] : 

AGl~ 9r = 5550 K - 2.31 T (87) 

rather than that used in [29J, and we also varied not only v\ but also two next constants, v<i 
and t>3. The v n values obtained are presented in the last line of Table 1. Variations of our 
V2 and v% with respect to their "purely deformational" values and v f lie certainly within 
the real accuracy of the original Blanter model as, first, this model disregards "chemical" 
contributions to i>2 and v% which can be quite notable (which is illustrated, in particular, 
by comparison of results of calculations of C-C interactions in ferrite based on ab-initio 
[20] and purely deformational [37] approaches [33]) and, second, it neglects both possible 
contributions of not-nearest Kanzaki forces [29j and a probable variation of phonon spectra 
with temperature (which was not measured in 7-iron but is very pronounced in the BCC 
iron [34, 35J). In Figs. 1 and 2 we present the carbon activity a^f gr (xq,T) and the 
equilibrium phase diagram ferrite-austenite calculated using the PCA expression (I77|) for 
A with our v n from Table 1, together with experimental data and the results of calculations 
based on the phenomenological model by Agren [36J. 
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Figure 3: Illustration of our method of estimates of saddle-point interactions A n using 
interpolation of configurational interactions v n = v(R n ). Circles show values of v n , and 
triangles, values of An obtained as described in the text. 

Figure 4: Values of carbon-carbon interactions used in the present work. Open circles: 
configurational interactions v n = v(R n ); black triangles: saddle-point interactions A„ = 
A(R^). Dashed lines connect the neighboring A n values to guide the eye. 

Table 2. Saddle-point interactions A n = A(R* P ) (in kelvin) of carbon atoms in austenite 
for vectors R n p = (R n — R sp ) where R sp is the saddle-point position of carbon atom 
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K p /Ri 


0.87 
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1.12 
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1.32 
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1.5 
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1.66 
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1.80 
12 


1.94 

8 


2.06 
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2.18 
12 


2.29 
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A„ 


1470 


1336 


1228 


229 


-924 


-929 


-543 


133 


133 


286 
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11 


12 


13 


14 


15 


16 


17 


18 


19 




R-nl Rl 
7 sp 


2.40 

8 


2.50 
14 


2.60 
16 


2.70 
4 


2.78 
16 


2.87 
16 


2.96 

8 


3.04 
20 


3.12 

8 




A n 


622 


564 


144 


-160 


-310 


-269 


34 


-43 


-59 





Let us now discuss the saddle-point interactions A n = A(R* P ) where R* p = (R^— R sp ), 
and R sp is the saddle-point position of carbon atom. In the second and the third lines of 
Table 2 we show the first 19 distances R s ^ = |R* P | and the coordination numbers that 
correspond to these R n p . To illustrate the distribution of vectors R^ p in the FCC lattice, 
below we present the values of components of lattice vectors R^ = (R^ p + R sp ) (in a/2 
units) for the first eight coordination spheres of the point R sp =(0.5,0.5,0): 

R; = (0, 1, ±1), (1, 0, ±1); R 2 = (0, 2, 0), (2, 0, 0), (1, 1, 0), (1, 1, 0); 
R^ = (2, 1, ±1), (1, 2, ±1), (0, 1, ±1), (1, 0, ±1); 

R^ = (2, 2, 0), (1, 1, 0), (0, 0, ±2), (1, 1, ±2); R' 5 = (2, 1, ±1), (1, 2, ±1); 

R^ = (0, 2, 0), (2, 0, 0), (1, 3, 0), (3, 1, 0), (0, 2, ±2), (2, 0, ±2), (1, 1, ±2), (1, 1, ±2); (1, 1, ±3); 

R' 7 = (0, 3, ±1), (3, 0, ±1), (2, 1, ±1), (1, 2, ±1), (3, 0, ±1); 

R 8 = (2, 2, 0), (2, 2, 0), (3, 1, 0), (1, 3, 0), (0, 1, ±3), (1, 1, ±1). (88) 

Models for estimation of the saddle-point interactions A n can be constructed similarly 
to those for the configurational interactions v n . The long-ranged deformational contribu- 
tions to A n can be calculated using the general expression for deformational interactions 
with the arbitrary Kanzaki forces discussed in [29] . while the short-range chemical con- 
tributions can be estimated treating several first A n as adjustable parameters, as it was 
made for the v n . However, in this work we restrict ourselves by illustrative estimates of 
A n based on some interpolations between v n values in Table 1 and several simple assump- 
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tions. First, we assume that both the chemical and the deformational contributions to 
A n depend only on the distance = |R^ P | and vary with R s v smoothly Second, we 
assume that for short distances < i?2=l-41i?i, the A n values are mainly determined 
by the chemical contributions, while for longer Rff > R 2 , these values are mainly de- 
termined by the deformational contributions. Third, we assume that the dependence of 
the configurational interactions v n on R n is similar to the dependence of the saddle-point 
interactions A n on for both chemical and deformational contributions. Then the 
"chemical" interactions Ai, A2 and A3 can be estimated using the linear interpolation 
between v\ and v 2 values, as shown in Figs. 3 and 4. For the "deformational" A n with 
n > 4 or R^ > R 2 , the analogous estimate of A n includes the following two steps: 

(A) Interpolation of dependence v(R) using v n values in Table 1, which yields "pre- 
liminary" values shown in Fig. 3. 

(B) Scaling of these Ai°' by some factor a, 

A n = aAg\ (89) 

with the value a determined from the fit of the diffusivity D calculated according to Eqs. 
(J75]l - (j75|l to the experimental data about the diffusivity of carbon in austenite. 

The first step (A) is illustrated by Fig. 3. This figure shows that to obtain an adequate 
interpolation v(R), the regions of long and "intermediate" distances R should be treated 
differently. At long distances R > R4, we can use the simple linear interpolation between 
neighboring v n values, while at R 2 < R < R4, some smooth curve should be drawn 
between v 2 , t>3 and V4 values. For these intermediate R, we interpolated v(R) by a simple 
power law: 

R 2 <R<R 3 : v(R) = C 2 (R 3 - R) m , 

R 3 <R<R 4 : v(R) = C 4 (R - R 3 ) m (90) 

where constants C 2 and C4 are determined by the conditions: v(R 2 ) = v 2 , v{R$) = V4. 
For the exponent m we tried two values, 2 and 4, and the value m = 4 was found to be 
more suitable for the fit mentioned in the point (B). The resulting interpolation v(R) is 
shown in Fig. 3 by the dashed curve. 

For the second step (B), the physical arguments in favor of the model relation (I8"9~j) 
can be seen from the general expression for deformational interactions given by Eq. (11) 
in [29]. According to this expression, the deformational interaction V$ = V d (Ri — Rj) 
between two atoms positioned at Rj and R^ is proportional to the integral over wave- 
numbers k in the Brillouin zone of some expression that includes the product of two 
appropriate Kanzaki forces, f£ and f^, while each of these forces is proportional to the 
amplitude of displacements of neighboring host (iron) atoms due to the presence of an 
impurity (carbon) atom at site i or j. Therefore, for the configurational interactions v n , 
the deformational contributions are proportional to the product of two Kanzaki forces, 
f° p and f° p , that describe the displacements of iron atoms induced by a carbon atom 
positioned in the octo-pore. At the same time, for the saddle-point interaction A n , one of 
these factors is replaced by a Kanzaki force f* p that describes the analogous displacements 
of iron atoms but induced by a carbon atom in the saddle-point position, for which the 
carbon-iron distance RFe-c is by l/y/2 times smaller than that for a carbon atom in 
an octo-pore. Therefore, this Kanzaki force ff p can be expected to notably exceed f° p . 
Hence the factor a in ( I89p . which qualitatively describes the relative scale of deformational 
contributions to the v n and to the A n values, can notably exceed unity. 
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Figure 5: The diffusivity Dc(xc,T) of carbon in austenite. The symbols denote exper- 
imental values from [38] presented in [3]. The solid lines are calculated using Eq. fP75|) 
and the values v n and A n from tables 1 and 2. The dashed lines show the results of 
calculations by Agren [5] based on his phenomenological model. 

The description of available experimental data about the diffusivity Dc(xq,T) [2E] 
by our model with the choice a = 2.9 in Eq. (18*91 is shown in Fig. 5. This description 
corresponds to the following values of the saddle-point energy E ac , the product ujq a 2 , 
and the frequency u>q* in Eqs. (!75|) and (!76j) : 

E ac = 17700 K, (91) 
uj e c ff a 2 = 0.225 cm 2 /sec, u e c ff = 1.76 • 10 14 sec" 1 (92) 

where the value a = 3.58 A for 7-iron [ID] is used. Let us note that the choice a = 3.04 in 
Eq. (189]) would yield the values of Dc{xc, T) which virtually coincide with those obtained 
by Agren [3]. However, the choice a = 2.9 seems to better describe the low-temperature 
data by Wells et al. [38J which agree with those obtained by Smith [39]. The saddle-point 
interactions A n which correspond to a = 2.9 are presented in Table 2 and Fig. 4. 

Both experimental and theoretical results presented in Fig. 5 show that the diffusivity 
Dq sharply rises with increase of the carbon concentration xc- In accordance with Eqs. 
fTTOj) . (TTTT) . (PT51) and fPTBl) . it seems to indicate on the presence of a significant attraction 
in the saddle-point interactions A n , as this attraction lowers the saddle-point energy Eq P 
for the inter-site jumps of carbon atoms. In the model estimates of interactions shown in 
Fig.4, it corresponds to the presence of significant negative A n at "intermediate" carbon- 
carbon distances R lying in the interval 1.6i?i < R < 2R±. The rise Dc(c,T) with xc 
mentioned above can imply that such significant attraction is present not only in our model 
estimates, but also in the real saddle-point interactions of carbon atoms in austenite. 

Let us now discuss the values of the pre-factor cu^ and the "transition state entropy" 
ASq P in ([S]) which correspond to the estimate (I9"2"l) . The attempt frequency co pv in §E§ for 
our case can be estimated as the frequency uJq of local vibrations of carbon in austenite. 
These vibrations have been experimentally studied by Sumin et al. [41J who found: 

col - 75meV = 1.14- lO^sec^ 1 . (93) 

Note that this Uq exceeds the Debye frequency of 7-iron, uP D = 0.43 • 10 14 sec -1 [12], by 
about three times. Then using Eq. ([8]) with u pv =Uq and uj^/=lOq from (1921) . we obtain: 

ASg p ~0.4, cu s c p -0.9 ul. (94) 

These relations show that the "softening" of saddle-point frequencies uj S q with respect 
to uPq for carbon in austenite is rather weak (if any), unlike Fe-Cu substitution alloys 
discussed in Sec. 2.1, while the saddle-point entropy ASq P is by an order of magnitude 
lower than the analogous AS^ P V and AS§^ values for the Fe-Cu alloys (as estimated 
by SF [2]). The difference can be related (at least, partly) to the above-mentioned 
inequality u^, ^> iS^ which implies that the dynamics of carbon atoms in austenite is 
much faster than the iron atom dynamics. In such conditions, an assumption of "a local 
thermodynamic equilibrium" for the saddle-point transition state, as well as the entropy 
notion for this state can be not fully adequate and should be used with some caution. 

Note that Eqs. (1941) correspond to the pre-factor of the factor 7 in (1761) which 
determines the diffusivity (1751) in the dilute alloy limit. Therefore, these equations have no 
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relevance to the illustrative estimates of carbon-carbon interactions discussed above, but 
they provide some definite information about the microscopic characteristics of diffusion 
of carbon in 7-iron. 

Finally, let us compare the microscopic description of thermodynamic and diffusional 
characteristics of carbon in austenite presented in this work to their phenomenological 
description developed by Agren [2] E] EE]- Both approaches use a similar number of 
adjustable parameters, and the quality of agreement between the results obtained and 
the experimental data shown in Figs. 1, 2 and 5 is similar. However, the microscopic 
approach seems to provide a more adequate physical understanding of the phenomena 
considered. It also opens possibilities for developments of fully first-principle descriptions 
with no adjustable parameters, as demonstrated by SF [TJ] for Fe-Cu alloys. In addition 
to that, the microscopic expression (175]) for the diffusivity seems to elucidate a number of 
principal points not discussed earlier. First, it shows that the diffusivity can be written 
in the form of the product of "thermodynamic" and "kinetic" (or " saddle-point" ) factors, 
and the thermodynamic factors include not only the so-called Darken-factor d In ax/d In c 
usually written [2], but the concentration derivative of the activity itself, dax/dc. Second, 
the microscopic relations (T75]) -f r78]) enable us to estimate the "transition state entropy" 
AS SP from experimental data, as was demonstrated for carbon in austenite. Third, these 
microscopic relations allow to relate the concentration dependence of the activity ax and 
the diffusivity Dx to both the configurational and the saddle-point interactions between 
interstitial atoms X, in particular, between carbon atoms in austenite. Therefore, the 
analysis of experimental data about ac(xc,T) and D c (xq,T) enables us to get an idea 
about the type and the scale of these interactions. 

5 CONCLUSIONS 

Let us summarize the main results of this work. The fundamental master equation for the 
probability of various atomic distrubutions in an alloy is used to derive the basic equations 
of diffusional kinetics in alloys. The microscopic parameters entering these equations can 
be calculated by ab-initio methods, as was demonstrated by SF for iron-copper alloys 
[2] , or using various theoretical models. For substitution alloys, the diffusional transfor- 
mation kinetics is described by the "quasi-equilibrium" kinetic equation (QKE) derived 
in Sec. 2.1. This equation ([28]) generalizes its earlier version presented in [15] by taking 
into account possible "interaction renormalization" effects which can be important for the 
vacancy-mediated kinetics [TT1[T3]. In sec. 2.2 we describe the calculations of local chemi- 
cal potentials Aj and correlators bij entering the QKE (125]) with the use of some analytical 
methods which combine simplicity of calculations with a high accuracy, particularly for 
dilute alloys. In Sec. 2.3 we reduce the QKE (128]) describing the vacancy-mediated ki- 
netics to the kinetic equation for some equivalent direct-atomic-exchange model which is 
suitable for computer simulations. 

The microscopic equations describing diffusion of interstitial atoms X in an interstitial 
alloy Me-X are derived in Sec. 3. These equations have a simple form (I65p or (168 p . which 
enable us to obtain the explicit analytical expressions for the diffusivity D = Dx- These 
expressions for D given by Eqs. (T7T]) or (175]) have a simple form of products of three 
factors: the concentration derivative of the thermodynamic activity ax of atoms X; the 
correlator bx given by Eqs. (IBTj) or (175]) which describes the influence of interactions 
between atoms X on the activation barrier for the inter-site jumps of atoms X; and the 
concentration- independent factor 7 describing the diffusivity in the dilute alloy limit. This 
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microscopic expression for D differs notably from those used in phenomenological treat- 
ments PQ-[I], in particular, by the presence of the concentration derivative dax/dc rather 
than the so-called "Darken-factor" d In ax/d In c usually written. We also derive equations 
describing diffusion of interstitial atoms X in a multicomponent alloy (MeiMe2...)-X. 

In Sec. 4 we apply the results of Sec. 3 to microscopically treat the problem of diffu- 
sion of carbon in austenite discussed by a number of authors [H [21 El II]- Our treatment 
is based on the microscopic model of C-C interactions in austenite suggested by Blanter 
[16] which supposes a strong "chemical" repulsion at short C-C distances R and a purely 
deformational interaction at longer R. To estimate the configurational interactions v (R) 
which determine the carbon activity ac, and the "saddle-point" interactions A(i?) which 
determine the above-mentioned correlator bx=bc, we use some plausible assumptions 
about the dependences v(R) and A(i?) which include adjustable parameters. The inter- 
action models obtained enable us to describe both the thermodynamic and the diffusional 
properties of carbon in austenite at the same level of accuracy as that achieved in phe- 
nomenological treatments [TJ [21 El IH EE]- At the same time, the microscopic approach 
used enables us to make a number of qualitative conclusions about the carbon-carbon in- 
teractions and the characteristics of diffusion of carbon in austenite, in particular, about 
the presence of a significant C-C attraction at intermediate R and about a rather low 
value of the "transition state entropy" ASq P given by estimate (194P . 

The authors are much indebted to F. Soisson for numerous valuable discussions, as 
well as to G. F. Syrykh, for drawing our attention to the important experiments [JT]. The 
work was supported by the Russian Fund of Basic Research (grant No. 09-02-00563); by 
the fund for support of leading scientific schools of Russia (grant No. NS-7235.2010.2); 
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